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Context. GS 1826-238 is a well-studied X-ray bursting neutron star in a low mass binary system. Thermal Comptonisation by a hot 
electron cloud (kT e ~ 20 keV) is a widely accepted mechanism accounting for its high energy emission, while the nature of most 
of its soft X-ray output is not completely understood. A further low energy component is typically needed to model the observed 
spectra: pure blackbody and Comptonisation-modified blackbody radiation by a lower temperature (a few keV) electron plasma were 
suggested to explain the low energy data. 

Aims. In order to better characterise the nature of the low energy emission and the bolometric output of the source, the steady 
emission of GS 1826-238 is studied by means of sensitive, broad band (X to soft Gamma-rays) measurements obtained by the 
INTEGRAL observatory. 

Methods. In this data analysis, the newly developed, up-to-date Comptonisation model comptb is applied for the first time to study 
' effectively the low-hard state variability of a low-luminosity neutron star in a low-mass X-ray binary system. Three joint ISGRI/JEM- 

X data sets (two from observations performed in 2003 and one from 2006) were analysed. 

Results. We confirm that the 3-200 keV emission of GS 1826-238 is characterised by Comptonisation of soft seed photons by a 
hot electron plasma. A single spectral component is sufficient to model the observed spectra. At lower energies, no direct blackbody 
emission is observed and there is no need to postulate a low temperature Compton region. Compared to the 2003 measurements, the 
plasma temperature decreased from ~ 20 to ~ 14 keV in 2006, together with the seed photons temperature. The source intensity was 
also found to be ~ 30% lower in 2006, whilst the average recurrence frequency of the X-ray bursts significantly increased. Possible 
explanations for this apparent deviation from the typical limit-cycle behaviour of this burster are discussed. 
Conclusions. 
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fl") 1. Introduction in the burst decay was also reported (Thompson et al., 2005). 

I " The lack of evidence for photospheric radius expansion during 

> ■ The persistent X-ray source GS 1826-238 was discov- ±e burstfj (jn >t Zand et ^ lg9g) and j considerations on 

^ dUnng , a 1 ^o a o n ° e T UV ? f ft G T, Satelhte m , 19 u 88 the optical counterpart (LMXB, short orbital period) limit the 

K£ , (Makino et al., 1988). Early band-limited observations of this yalue of ±e distance tQ ±e SQWX in the ^ 4 _ 9 k ran 

^■object were performed by TTM (in t Zand, 1992), and (Barret et al 2000) 
' ROSAT (Barret, Motch & Pietsch, 1995) at low energies and 
by CGRO-OSSE (Strickman et al., 1996), in the hard X- The hlstor y of sensitive, wide band measurements of 
rays. An optical counterpart in a low inclination system the stead y emission of GS 1826-238 includes results from 
and with a likely orbital period of 2.1 h was reported BeppoSAX (in 't Zand et al., 1999, Del Sordo et al., 1999), 
(Homer, Charles & O'Donoughue, 1998). The detection of type RXTE (Barret et al., 2000), RXTE plus Chandra 
I X-ray bursts (Ubertini et al., 1997) by BeppoSAX ultimately (Thompson et al., 2005), and RXTE plus XMM 
established that the source is a neutron star (NS) in a low mass (Thompson et al., 2008). Since the BeppoSAX era, the source is 
binary system (LMXB). X-ray bursts were found to recur every alwa y s detected in a hard (low) spectral state, its high energy 
few hours in a quasi-periodical fashion (Ubertini et al., 1999). emission extending well in excess of 100 keV. A simple, 
What makes this burster unique among its class is the remark- empirical, exponentially-cutoff power law model with E c ~ 
able stability of the recurrence over a time-scale of several years 50 keV y ields acceptable fits to the high energy data (e.g. 
(Cocchi et al., 2001). Bursting behaviour is clearly related to the Del Sordo et 19 ") but whenever more physical models 
persistent emission (Cornelisse et al., 2003), the event frequency are considered, the spectra are generally well described by 
increasing with the measured X-ray flux. RXTE data permitted a Comptonisation scenario, where soft seed photons are up- 
more accurate studies of the burst recurrence properties (see scattered to high energies by a hot electron plasma (kT e ~ 20 
e.g. Galloway et al., 2004), and a possible 611 Hz oscillation keV )- Alon 8 wlth Comptonisation, either a low energy black- 

body (BB) spectral component (in 't Zand et al., 1999) or a 

Send offprint requests to: M. Cocchi low-temperature Compton component (Thompson et al., 2005) 

* e-mail: Massimo.Cocchi@iasf-roma.inaf.it are included to fit the data. The detection of a ~ 6 keV iron line 
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Table 1. Observation log. 



INTEGRAL revolution 


061-064 


119-122 


495 


obs. year 


ZUUj 


zUUj 


zuuo 


obs. day 


Apr 15-23 


Oct 06-15 


Nov 2 


JEM-X live time (ks) 


21.9 


32.8 


49.9 


ISGRI live time (ks) 


134.4 


206.5 


50.3 


# of bursts (JEM-X) 


3 


5 


4 


average burst delay (h) 


3.79 


3.77 


3.40 


ASM intensity (mCrab) 


37 + 7 


40 + 8 


25 ±4 



" average during the INTEGRAL observations, 2-10 keV. 



was reported by Barret et al., 2000, who also included reflection 
to fit their data. Evidence for transient soft blackbody emission, 
possibly related to irregularities in the time behaviour of the 
X-ray bursts, has recently reported by Thompson et al. (2008). 

We present wide band spectral analysis of GS 1826-238 per- 
formed by INTEGRAL data. In the next section, observation and 
data analysis details will be given. In Sect. 3 we present and dis- 
cuss our results along with our best fit models, while the implica- 
tions of these results and their relation to the current knowledge 
of this peculiar X-ray source are debated in Sect. 4. 

2. Observations and data analysis 

The public database 1 of INTEGRAL (Winkler et al., 2003) in- 
cludes a big set of useful pointings of GS 1826-238, thanks 
to the observing policy of the telescope which includes peri- 
odic scans of the Galactic plane and the Galactic bulge region 
(GS 1826-238 is only ~ 10 degrees from SgrA*). 

Detailed analysis of the whole 

INTEGRAL GS 1826-238 public database, focused on the 
source long term monitoring, is to be presented elsewhere 
(Cocchi et al., in prep.). The present work, focused on wide 
band spectroscopy of the persistent emission, relies on data 
obtained by the imaging detectors JEM-X (Lund et al., 2003) 
and IBIS/ISGRI (Lebrun et al., 2003, hereafter ISGRI), which 
operate in the 2-30 keV and 17-300 keV passbands, respec- 
tively. 

Being the source variable, though always observed in a 
low-hard spectral state, homogeneous data sets are needed: 
this means that one has to choose simultaneuos JEM-X+ISGRI 
pointings of sufficient length to get good statistics and as close 
in time as possible. Due to the different field of view of the 
two detectors and the typical dithering observing pattern of 
INTEGRAL, JEM-X live time is considerably less than the 
ISGRI one. 

Taking into account the above limitations, we selected 
three data sets, the first one belonging to the satellite revo- 
lutions 061-064 (visibility window of Spring 2003), the sec- 
ond one to the revolutions 119-122 (Autumn 2003), and the 
third one to revolution 495 (Autumn 2006). Together, the three 
sets account for 104.6 ks and 391.2 ks live time (JEM-X 
and ISGRI, respectively). Analysis of RXTE/ASM monitoring 
data of GS 1826-238 obtained in the same epochs of our 
INTEGRAL observations indicate that the 2-10 keV source in- 
tensity significantly varied from 2003 to 2006, being - 40 
mCrab in 2003 and only ~ 25 mCrab in Autumn 2006 (see Table 
1). This is confirmed with higher accuracy and on a wider energy 
band by the present measurements (see Sect. 3). 

1 http://isdc.unige.ch/index. cgi?Data+info 



In order to avoid possible low-energy spectral contamination 
from the X-ray burst emission, the time intervals of the observed 
bursts, each event typically lasting ~ 100 s in JEM-X data, were 
removed from the good time of the data sets. Details about the 
analysed data sets such as live time, number of bursts detected 
and inferred average burst periodicities, are summarised in Table 
1. 

JEM-X data were reduced and extracted by the stan- 
dard INTEGRAL off-line analysis (OSA) pipeline 2 release 7.0. 
Presence of systematic effects at the lowest energies suggested 
a discarding of the spectral bins in the 2-3 keV range. ISGRI 
data were reduced by the standard OSA 6.0 pipeline, while spec- 
tral extraction was performed by a dedicated software based on 
cross correlation, which is freely available to the community 3 . 
The effectiveness and the limitations of the software in terms 
of calibration and systematics have already been discussed else- 
where (Segreto & Ferrigno, 2006; Mineo et al., 2006). Finally, 
the joint JEM-X/ISGRI spectra were analysed and best-fitted by 
the standard software tools of the Xspec package 4 , version 12. 

Due to cross calibration uncertainties between the two instru- 
ments, whose effects are expected to vary with the observation 
epoch 2 , a constant normalisation factor was included in each fit. 
We kept ISGRI as a reference, allowing the JEM-X constant to 
vary. No systematics, accounting for possible instrumental mis- 
calibration, had to be included for both JEM-X and ISGRI. 

Our data cannot constrain the low energy absorption along 
the line of sight, so the Nn value, modeled by Wisconsin ab- 
sorber in Xspec , was kept fixed to 3 x 10 21 H cm 2 , which 
is the best determined value from recent joint XMM-RXTE 
measurements (Thompson et al., 2008). Slightly different val- 
ues, 5 x 10 21 H crrT 2 and 1.5 x 10 21 H ctrT 2 were determined 
respectively by previous ROSAT and BeppoSAX measurements 
(Barret, Motch & Pietsch, 1995; in 't Zand et al., 1999). All the 
values are of the same order of magnitude of the interpolated one 
(Dickey & Lockman, 1990) that can be derived by nh, a widely 
used ftool 5 . We already verified that these slight differences in 
the /Vh value do not affect the parameters of our best fits and the 
possible undetection of further soft components (see Sect. 3). 

3. Spectroscopic results 

Our wide band INTEGRAL spectra of GS 1 826-238 can be satis- 
factorily fit in a first approximation by a simple, empirical, expo- 
nentially cutoff power-law model, Xspec 's cutoffpl (see Table 
2). The power-law photon index is ~ 1 .6 and the cutoff energy is 
~ 50 keV in all the observations. Such values are in good agree- 
ment with all the previous results based on this simple model, 
thus suggesting that the overall spectral shape (and so the spec- 
tral state) of the source is almost constant with time on a time 
scale of several years. 

From a less empirical, more physical point of view, a single 
high energy cutoff power law is a clear indication of a single 
Compton region at work. This is (at least partially) in agreement 
with most of the previous results, but not with those obtained 
by Thompson et al. (2005) (see also Thompson et al., 2008), 
who included a second, distinct Compton component with a 
different geometry and electron temperature to fit their joint 
Chandra /RXTE data: in our case, as cutoffpl is an acceptable 
model, there is a clear indication that the data themselves do not 

2 http://isdc.unige.ch/index.cgi?Soft+info 

3 http://www.ifc.inaf.it/~ferrigno/INTEGRALsoftware.html 

4 http://heasarc.nasa.gov/docs/xanadu/xspec/index.html 

5 http://heasarc.nasa.gov/docs/software/ftools/ftools_menu.html 
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Fig. 1. The INTEGRAL spectrum of GS 1826-238 as observed in Autumn 2003 (brighter spectrum) and Autumn 2006. Simple crosses and dotted 
crosses are used for ISGRI and JEM-X, respectively. Residuals to the best comptb fits are shown in the lower panels for the 2003 (above) and 2006 
data (below). 



require the use of too many free parameters to explain the ob- 
served emission from GS 1826-238 . 

Taking all this into account, we tried to fit the data with 
a single thermal Compton component. The recently developed 
thermal plus bulk Comptonisation model (Farinelli et al., 2008), 
comptb , included in the standard Xspec distribution 6 , is 
presently the most up-to-date Comptonisation model in a dif- 
fusion regime (such approximation holds well for NS-LMXBs). 
comptb is a general model including, in a self-consistent way, the 
seed photons emission together with the Comptonised emission, 
whether the latter is either purely thermal (as in the well-known 
model comptt ), or mixed dynamical (bulk Comptonisation) 
plus thermal. As a detailed description of the model is given in 
Farinelli et al. (2008), we here only point out that the free param- 
eters of comptb are the blackbody seed photon and the electron 
temperatures (kTs and kTe, respectively), the Comptonisation 
index a, the bulk parameter 6, the covering factor log(A) and 
the normalisation. 

As already mentioned above, the spectrum of 
GS 1826-238 is well described by a simple function (cut- 
offpl), which mimics unsaturated Comptonisation for static 
plasma and no blackbody-like peak is observed at low energies: 
the case of the seed blackbody radiation dominated by the 
Compton emission is modeled by comptb with a high value of 
the covering factor, i.e. A » 1. Moreover, GS 1826-238 be- 
ing a low luminosity (a few % Eddington) NS-LMXB in a 
low (hard) state, the detection of a bulk emission is unlikely 

6 http://heasarc.nasa.gov/docs/xanadu/xspec/newmodels.html 



(Paizis et al., 2006, see also Sect. 4.1). These two particular 
cases are well reproduced in our fits: when allowed to vary, 
both the log(A) and the 6 parameters are eventually pegged 
by the fitting procedure to their hard limits corresponding to 
total coverage (log(A)= 8) and no bulk emission (S = 0), 
thus reproducing a simple thermal Comptonisation spectrum 
without direct emergence of the underlying seed photons. 
We then fixed 5 and log(A) to the above limit values in all 
the fits, leaving only kTs, kTe and a as variable parameters. 
In this particular case, the dimension of the free-parameter 
space is identical to that of the well consolidated and widely 
adopted comptt model (Titarchuk, 1994), and the two models 
are substantially equivalent. 

The results of our comptb fits are reported in Table 2 for each 
of the three different observation periods. Spectra and residuals 
to the best fits are also shown in Fig. 1 for two selected obser- 
vations. A single comptb component with kT e in the range 15- 
20 keV effectively describes the shape of the spectra. Though 
the source persisted in its low-hard state, significantly different 
best fit values were obtained for both electron plasma and soft 
seed photon temperatures in the 2003 and 2006 observations. 
The 2006 data set required lower temperatures and a lower value 
of the index a, indicating an increased efficiency of the Compton 
upscattering. The 2-200 keV intensity also dropped by ~ 30% 
with respect to 2003. We carefully verified that the different val- 
ues of the main spectral parameters in 2006 are not biased by a 
probable long term variation of the JEM-X/ISGRI normalisation 
constant: in fact consistent kT e values are obtained even when 
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Table 2. Spectroscopic results. 



.0 F 



data set (revolution) 


061-064 


119-122 


495 






wabs*CompTB 


a 






kT s (keV) 
a index 
kT e (keV) 


1 06 +007 

1-05 !°1 
20.1 tg 


0.95 ± 0.05 
1.03 ±0.03 
19.4 +1 ° 


045 -III 

0.85 S 
13.5 ±0.8 


7 
£ 


const. b 

red. x 2 (d.o.f.) 


671 +0028 

"• u/ 1 -0.032 

1.076 (360) 


n 709 +0.027 
"•'J* _o.025 

1.092 (361) 


f) 0^9 +0-029 
-0.037 

1.094(344) 


> 


flux 2-20 keV c 


1.11 x 10-" 


1.34 x 10- 9 


1.02 x 10- 9 




flux 20-200 keV 


1.23 x 10-" 


1.34 x 10-" 


0.80 x 10-" 




flux 2-200 keV 


2.34 x 10~ 9 


2.68 x 10~ 9 


1.82 x 10~ 9 


l_J 


luminosity d (erg s -1 ) 


1.76 x 10 37 


1.92 x 10 37 


1.42 x 10 37 




r e 


3.4 


3.5 


5.3 




CAF d 


3.8 


3.9 


11.5 






wabs*CompTT 






kJ Q (keV) 


1.02 + 0.06 


0.91 ± 0.05 


0.33 ± 0.04 





t/ 



(keV) 



const. 

red. ^ 2 (d.o.f.) 
flux 2-200 keV 



21.5!" 

1, gl +0.33 

0.680 S 
1.064 (360) 

2.34 x 10-" 



20.0 !| 2 

4 IT +0.23 
-0.24 

756 +0026 
1.081 (361) 
2.68 x 10" 9 



14.6!!;? 

5.55 ! 033 

900 
u.you _ 004] 

1.070 (344) 

1.85 x 10" 9 



Ecu, (keV) 

r 

const. 

red. x 2 (d.o.f.) 
flux 2-200 keV 



WABS^CUTOFFPL 

51.9 + _ts 



1.57 ± 0.05 
0.784 !X 

1.022 (361) 
2.46 x 10" 9 



55.9 _2 g 



1.62 ± 0.03 

829 +0 028 

u.oz.y _ Q27 

1.065 (362) 
2.83 x 10" 9 



49.0 + f J 
1.62 ± 0.04 
0-825 $™ 
1.001 (345) 
1.91 x 10" 9 



" name of the Xspec model. n H fixed to 3.0 x 10 21 H cm -2 in all fits. 
All errors are 90% probability. 
b JEM-X/ISGRI normalisation ratio (see Sect. 2). 
c all fluxes in erg cirr 2 s~' . 

d unabsorbed, 0.1-200 keV, assumed distance of 7 kpc 
derived parameter (see Sect. 3) 



■ f assumed spherical geometry 



forcing the 2003 cross-calibration in the 2006 data set, and vice 
versa. Moreover, considering only ISGRI data for the fits does 
not alter the main parameter values. Finally, the RXTE /ASM 
data independently confirm in a narrower passband the drop in 
the flux (see Table 1). 

Best-fit excesses above ~ 75 keV are observed in the 2006 
data set (see lower panel of Fig. 1) and partly in the Autumn 
2003 one. But they are not detected in Spring 2003 and moreover 
their significance is poor. So we regard such residuals as random 
fluctuations. 

Interesting additional physical information can be extracted 
from the best-fit parameters of comptb ; for example, the 
cloud optical depth r can be derived by the eqs. [17]-[24] in 
Titarchuk & Lyubarskij (1995). The r values, assuming a spher- 
ical geometry, are included in Table 2. For a direct compari- 
son, the results obtained replacing comptb with the more popular 
model comptt (spherical geometry assumed) are also reported in 
Table 2. As expected, the main comptt parameter values are in 
good agreement with the comptb ones. 

Moreover, the order of magnitude of the radius of the (as- 
sumed) spherical region emitting the blackbody seed photons 
can be estimated: as the normalisation in comptb is defined as 
Cn = L39 /I>i , where L39 is the seed photon spectrum luminos- 
ity (in units of 10 39 erg s -1 ) and D10 is the source distance (in 




Fig. 2. Comparison of wide band spectra of GS 1826-238 ob- 
tained by different observations: our INTEGRAL spectra, as 
in Figure 1 (simple lines); RXTE (Barret et al., 2000, dashed 
line); Chandra +RXTE (Thompson et al., 2005, dash-dotted); 
BeppoSAX (in 't Zand et al., 1999, dots). The Fe lines, even when 
detected, are not included. 



units of 10 kpc), one can derive applying the Stefan-Boltzmann 
law the emitting region radius as 



Rbn = 0.88 x 10 



2 CH 2 D W 



(1) 



We obtain for the three observational data sets the values R^ = 
6, 8 and 24 km, respectively (a source distance of 7 kpc is as- 
sumed). 

The X-ray energy amplification caused by the Compton 
scattering can also be determined by computing the Compton 
Amplification Factor (CAF), defined as the ratio of the 
Comptonised bolometric energy flux to that of the blackbody 
seed radiation. The obtained values of the parameter for each ob- 
servation of GS 1826-238 are also listed in Table 2. Increased 
Compton efficiency, as already testified by the correspondingly 
decreasing value of the a index, is observed in 2006. 

4. Discussion 

4.1. GS 1826-238 in a unifying framework 

We applied a recently developed model for Compton upscat- 
tering, comptb , to study the low-hard state variability of 
GS 1826-238 . As Paizis et al. (2006) demonstrated, general 
Comptonisation models including bulk emission are very effec- 
tive in modeling the variety of NS-LMXBs spectral states and 
the state transitions in individual sources. This is obtained, in 
a natural way, by smooth changes in the relative importance of 
bulk over thermal Comptonisation. In this scenario having the 
interplay of bulk and thermal Compton emission as unifying pa- 
rameter, a low-hard state is expected to originate at low accre- 
tion rates, where the coronal Compton cooling is far from be- 
ing efficient, thus leaving out a hot screen which extends radi- 
ally and spherically, thus shielding the NS surface and allow- 
ing the direct accretion disk emission to be cooler. Bulk motion 
Comptonisation can still occur in this case, as the NS gravity is 
not balanced by the radiation pressure, but the effect is expected 
to be weak and probably washed in the hot screen with negligible 
effects for a distant observer. 
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The present results confirm the above picture for the case 
of GS 1826-238 , as pure thermal upscattering completely ac- 
counts for the source emission. The relative extension of the 
Compton cloud determines whether it is possible to observe a 
fraction of the seed photons' population originating from the NS 
surface and/or from the inner accretion disk. In our spectra the 
soft blackbody seed photons are not directly detected, given the 
inferred log A » 1, and no further soft components are needed 
to fit the data. 

4.2. Wide band spectroscopy: extra spectral components 
and variability 

Undetection of additional soft components in the 
INTEGRAL data is in agreement with the results obtained 
by Strickman et al. (1996) and Barret et al. (2000) for the same 
source, but not with those of in 't Zand et al. (1999) and Del 
Sordo et al. (1999), who proposed an extra blackbody (kTbb of 
~ 3 keV and ~ 1 keV, respectively), in addition to a thermally 
Comptonised component (or a Compton-like, in the case of 
Del Sordo et al., 1999), to fit the BeppoSAX data. Our upper 
limits for a further blackbody-like component were calculated, 
assuming a typical colour temperature of 1 keV and a distance 
of 7 kpc. We obtained blackbody emitting regions with radii 
R< 1.0 km and R < 0.7 km for the 2003 and 2006 data sets, 
respectively (90% confidence). 

Applying a pexrav model to our data as in Barret et al. (2000) 
does not improve the fits. The simple cutoff power law is repro- 
duced, but the reflection parameter is pegged to zero. This is 
expected, as cutoffpl always yields acceptable fits and there is 
no need for further hypotheses. 

For GS 1826-238 , Thompson et al. (2005) and (2008) pro- 
posed a model with two distinct Comptonising regions, one of 
spherical shape and a high (~ 20 keV) temperature, the other 
with disk geometry and a lower (~ 5 keV) temperature. The 
two separate regions had different opacities and seed photons 
populations. This two-cloud model actually fits our data, with 
values very close to those obtained by the above authors, but 
this is not very surprising because it has many free parame- 
ters. As a matter of fact, our data do not need too many free 
parameters because a simple cutoff powerlaw provides a satis- 
factory fit. Moreover, a single powerlaw is an indication of a 
single Comptonising region at work. In Thompson et al., 2005, 
much simpler models (e.g. cutoffpl) yielded acceptable fits to 
their data as well; interestingly, all the fits had very low (~ 0.7- 
0.8) reduced chisquared, probably indicating systematic effects 
at work, such as error overestimation or missed rebinning. 

It must be stressed that the presence of a second, softer, spec- 
tral component cannot be completely ruled out by the present 
INTEGRAL data. But its actual need should be justified by 
the overall X-ray spectral shape obtained including the lowest 
energy range (E< 2-3 keV), which is not covered by JEM- 
X. This is likely, as both the BeppoSAX and Chandra already 
published spectra start well below 1 keV. Nevertheless, both 
the previously reported ~ 3 keV blackbody-like component 
(in 't Zand et al., 1999) and the ~ 5 keV Comptonised emission 
(Thompson et al., 2005) are well inside the INTEGRAL pass- 
band and should be unambiguously constrained. This is not the 
case. 

The present analysis shows clear evidence for intensity and 
spectral variability of GS 1826-238 (even though almost unex- 
pected, given the peculiar long time-scale regularity of its burst- 
ing behaviour, which is suggestive of a stable, smooth accretion). 



In particular, significantly lower plasma temperatures and 2-200 
keV luminosity are observed in 2006 compared with 2003. 

A comparison of our results with previously reported wide 
band observations of GS 1826-238 is shown in Fig. 2. While 
at least approximately the canonical low-state spectral harden- 
ing with source dropping luminosity is well reproduced by our 
2003 data, the unusually low electron plasma temperature of 
2006 coupled with relatively low luminosity seems to contradict 
this consolidated scenario. Together with the probably erratic 
presence of blackbody-like emission at lower energies (see also 
below), the overall picture of GS 1826-238 is perhaps sugges- 
tive of the emergence of possible sub-states within its low-hard 
state. In such a spectral state, different observational characteris- 
tics could be attributed to differences in the configuration of the 
accreting system, with the size and geometry of the Compton 
cloud likely playing a pivotal role. 

4.3. X-ray bursts timing irregularities 

Taking into account the observed decrease of the X-ray bursts' 
wait times together with lower measured intensity (see Table 1), 
our November 2006 data show a clear deviation from the typical 
limit-cycle X-ray bursting behaviour of GS 1826-238 . Burst 
wait times AT are generally anticorrelated with very good ap- 
proximation to the source intensity Fx by the simple formula AT 
oc F^ 105 (see Galloway et al., 2004). If one were to consider the 
measured 2-200 keV flux as an indication of the actual accretion 
rate as well as the burst fuel accumulation rate, a deviation of ~ 
35% from the limit-cycle rule would be obtained with respect to 
e.g. the 2003 INTEGRAL data. Similar anomalies were already 
observed in the April 2003 RXTE data (Thompson et al., 2008). 
These results indicate that spectral measurements which do not 
include the very soft X-rays are unreliable tools to extrapo- 
late the bolometric luminosity, which can be widely underes- 
timated. In fact Thompson et al. (2008) reported of a simul- 
taneous XMM-observed transient blackbody-like spectral com- 
ponent. This very soft (a few tenths of keV) component ac- 
counted for the missing flux, allowing for a better estimate of 
the bolometric luminosity and so of the true accretion rate. Our 
INTEGRAL results can also be explained in the same natural 
way, leaving out early burst ignition and/or anisotropy consid- 
erations. But a further very soft component similar to the one re- 
ported by Thompson et al. (2008) would then remain undetected, 
being outside of the JEM-X passband. Nevertheless, our 2006 
data show a significant decrease of both the Compton cloud elec- 
tron and the seed photon temperatures, which is possibly related 
to an overall bolometric spectral variation of GS 1826-238 . 

The above reported erratic variability of the source, in terms 
of both transient soft components and/or Comptonisation pa- 
rameters, deserves to be studied in detail to be better under- 
stood and implies the need for further accurate wide band (from 
very soft X-rays up to soft Gamma-rays) monitoring programs 
of GS 1826-238 and other similar X-ray bursters. 
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